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ABSTRACT 

Recent studies suggest that Swift gamma-ray bursts (GRBs) may not trace 
an ordinary star formation history. Here we show that the GRB rate turns out to 
be consistent with the star formation history with an evolving stellar initial mass 
function (IMF). We first show that the latest Swift sample of GRBs reveals an 
increasing evolution in the GRB rate relative to the ordinary star formation rate 
at high redshifts. We then assume only massive stars with masses greater than 
the critical value to produce GRBs, and use an evolving stellar IMF suggested by 
Dave (2010) to fit the latest GRB redshift distribution. This evolving IMF would 
increase the relative number of massive stars, which could lead to more GRB 
explosions at high redshifts. We find that the evolving IMF can well reproduce 
the observed redshift distribution of Swift GRBs. 

Subject headings: cosmology: theory - gamma rays: bursts - stars: mass function 

1. Introduction 

Gamma-ray bursts (GRBs) are brief flashes of 7-rays occurring at an average detection 
rate of a few events per day at cosmological distances. Because of their very high luminosity, 
GRBs can be detected out to the edge of the visible Universe (Ciardi & Loeb 2000; Lamb 
& Reichart 2000; Bromm & Loeb 2002; Gou et al. 2004). Thus, GRBs are ideal tools for 
probing the star formation rate, the reionization history, and the metal enrichment history of 
the Universe (Totani 1997; Campana et al. 2007; Bromm & Loeb 2007). The advantages of 
GRBs over quasars for probing the high-redshift Universe had been discussed by Bromm & 
Loeb (2007). In addition, GRBs have been used as standard candles to constrain cosmological 
parameters and dark energy (Dai, Liang & Xu 2004; Friedman & Bloom 2005; Wang & Dai 
2006; Schaefer 2007, and references therein). 
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The association of long GRBs with core-collapse supernovae naturally suggests that the 
cosmic GRB rate should trace the star formation history. This gave rise to the expectation 
that GRBs may be a good tracer of cosmic star formation (Totani 1997; Wijers et al. 1998; 
Lamb & Reichart 2000; Blain & Natarajan 2000; Porciani & Madau 2001). However, it was 
found that the rate of GRBs increases with cosmic redshift faster than the ordinary star 
formation rate (SFR) does (Daigne et al. 2006; Le & Dermer 2007; Kistler et al. 2008, 2009; 
Yiiksel & Kistler 2007; Gen & Fang 2007; Li 2008; Wang & Dai 2009; Butler et al. 2010; 
Wanderman & Piran 2010). The reason for this discrepancy has been unknown. 

By investigating the redshift distribution of Swift GRBs, Guetta & Piran (2007) found 
that the observed high-redshift bursts are more than the expectation from an ordinary star 
formation history (SFH) and thus the high-redshift GRB rate is inconsistent with the one 
inferred from the current model for the SFR. Furthermore, Kistler et al. (2008) found that 
the GRB rate at redshift z ~ 4 is about four times larger than expected from star formation 
measurements. Daigne et al. (2006) concluded that GRB properties or progenitors must 
evolve with cosmic redshift to reconcile the observed GRB redshift distribution with the 
measured SFH. Li (2008) explained the observed discrepancy between the GRB rate history 
and the star formation rate history as being due to cosmic metallicity evolution, by assuming 
that long GRBs tend to occur in galaxies with low metallicities. However, very recently 
Levesque et al. (2010a,b) found several high-metaUicity long GRB host environments, which 
suggests that a low-metallicity cut-off is unhkely (also see Graham et al. 2009). Xu & 
Wei (2008) used a factitious stellar initial mass function (IMF) evolving with redshift to 
interpret the GRB redshift distribution. Gheng et al. (2010) suggested that this discrepancy 
could be eliminated if some high-redshift GRBs are ascribed to electromagnetic bursts of 
superconducting cosmic strings, although the existence of the superconducting cosmic strings 
has remained controversial. 

In this Letter, we first enlarge the GRB sample with 122 long GRBs observed by Swift. 
Then we interpret the latest Swift GRB redshift distribution using a reasonable evolving 
stellar initial mass function (IMF) proposed by Dave (2010). The structure of this paper is 
as follows: in section 2, we give an evolving initial mass function form, and in section 3, we 
show the analysis method. The results are presented in section 4 and conclusions are shown 
in section 5. 



2. An evolving initial meiss function 

The ordinary form of stellar IMF proposed by Salpeter (1955) is 

dN ^m-^-^^d\ogm. (1) 
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It was shown that the high-redshift GRB rate exceeds the expectation based on the above 
SFR. This leads us to invoking a different form of stellar IMF. The possibility of an evolv- 
ing IMF was discussed several times in the literature (e.g. Larson 1998, 2005; Ferguson, 
Dickinson & Papovich 2002; Fardal et al. 2007; van Dokkum 2008). Kroupa (2001) pointed 
out that a universal IMF is not expected theoretically, though no variations had been un- 
equivocally detected in the studies of local star-forming regions. Scalo (1998) mentioned 
that although the IMF index may vary at different redshifts in the Universe, its average 
value is close to the Salpeter value. Wilkins, Trentham & Hopkins (2008) independently 
determined the cosmic stellar mass growth rate by compiling observations of stellar mass 
densities from the literature and suggested an evolving IMF to interpret the discrepancy 
between the stellar mass density and SFH. In order to reconcile the discrepancy between the 
theory predicting the galaxy stellar mass-SFR relation with the observations, Dave (2008) 
proposed an evolving IMF with the following form 

dN f m~°-^ for m < rhiMF /o\ 

— = ^{m)QC-\ _i3 . ^ . (2) 

diogm [ m for m > miMF, 

where rhiMF = 0.5(1 -|- z^Mq, which has been constrained by requiring non-evolving star 
formation activity parameter. It is worth noting that this evolving IMF is only constrained 
out to 2; ~ 2 from the galaxy stellar mass-SFR relation, though its predictions are consistent 
with the other observations out to z ^ A. Chary (2008) found that this IMF would produce 
sufficient ionizing photons to account for late reionization of the intergalactic medium if it 
evolved out to z > 4. So the IMF suggested by Dave (2008) can be used to 2; > 4. More 
recently, Dave (2010) found miMF = 0.5(1 + z)^^^'"^^^ Mq using the Herschel data. 



3. The method 

We consider a spatial volume V at redshift z. The IMF can be written as A^lm), so 
J^' Am^{m)d\ogm = _Rsfr^) where -Rsfr is the SFR, mi is the largest mass of stars and 
is the smallest mass of stars. We consider only massive stars with masses larger than 3OM0 
can produce GRBs (Woosley 1993; Bissaldi et al. 2007)0, so 

-KgRB OC — - A — -jfni TTj -KSFR, [o) 



^The lower limit mass of a star that can collapse to GRB is uncertain at present. But this value is 
unimportant in our analysis below. The best fitting parameters will shift slightly when the lower limit mass 
is changed. But this evolving IMF could still interpret the GRB redshift distribution. 
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where K is a constant to be constrained and Rqrb is the rate of GRBs, representing the 
number of GRBs per unit time per unit volume at redshift z. We use the SFR derived by 
Hopkins & Beacom (2006), 



log Rsfr{z) = a + b log(l + z) , (4) 



with 



(-1.70,3.30), ^< 0.993 
(a, 6) = <( (-0.727, 0.0549) , 0.993 <z< 3.80 (5) 
(2.35, -4.46) , z > 3.80 

Le Borgne et al. (2009) used mid- and far-infrared observations to constrain the SFR and 
found that SFH is well-constrained and consistent with direct measurements from Hopkins 
& Beacom (2006). So we also use the results of Hopkins & Beacom (2006). 

Then the observed rate of GRBs within z z + dz and L L + dL is 

^ = m^^^^^dLdz, (6) 
dt ^ 'l + z Att dz ' ^ ^ 

where ^{L) is the beaming-convolved luminosity function of GRBs, {1 + z)~^ is due to 
cosmological time dilation and AQg = 1.4 sr is the solid angle covered on the sky by Swift 
(Salvaterra & Chincarini 2007). In a flat universe, the comoving volume is calculated by 



= 4^^com^ , (7) 

where the comoving distance is 

In the calculations, we use fi^ = 0.3, = 0.7 and Hq=70 km Mpc~^. There are many 
luminosity function forms in the literature. We use the Schechter-function form 

$(L) = i-(^|-yexp(-L/LO, (9) 

where /3 and are constant parameters to be determined by the observational data. The 
observed distribution of Ljso is then given by 



AN{L) = $(L) 



I 



— ; : ; — dz 

1 + z 47r dz 



ALAtobs, (10) 
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where 2;max = -2max(-Z^iso) IS the maximum redshift up to which a GRB with luminosity Liso 
can be detected by Swift, solved from equation Lii^{z) = L^so- 

The isotropic-equivalent luminosity of a GRB can be obtained by Liso = -E'iso(l + -2)/^9o 
(Kistler et al. 2008). We use 122 long GRBs observed by SwifQ till GRB 090726 (Butler 
et al. 2010). The distribution of Liso for the 122 GRBs in the sample is shown in Fig.l. 
The luminosity threshold can be approximated by a bolometric energy flux limit = 
1.2 X 10~*erg cm~^ s~^. The luminosity threshold is then 

Lum = ^riDlFii^ , (11) 

where -Dl is the luminosity distance to the burst. 

With the above luminosity threshold and an adopted GRB rate history, the observed 
luminosity distribution can be fitted by an intrinsic Schechter luminosity function with a 
power-law index /3 = —1.12, a characteristic luminosity = 9.16 x 10^^ erg s~^ and 
AtobsK = 69858.51 with Xr = 1-15. 



4. The results 

In order to study the rates of GRBs and star formation, it is convenient to use a 
dimensionless Q, where Q = Q{z) is defined by (Kistler et al. 2007) 

The coordinate Q is particularly useful in binning the data, as the definition of Q has taken 
into account both the effect of the comoving volume and the effect of cosmic time dilation. 
For example, when the comoving rate density of GRBs was a constant, in each equally sized 
bin of Q, the observed GRBs number would be a constant. The complete GRB selection 
function is very difficult to determine (Coward 2007). We choose the GRBs in the cuts 
Liso > 0.8 X lO^^erg s~^ and in the redshift range — 4 (Kistler et al. 2008). This method 
can reduce the selection effect by removing many low-^, low-Ljso bursts that could not have 
been seen at higher redshifts. There are 72 GRBs in this sample. The SFR fit from Hopkins 
& Beacom (2006) in this range is shown as the dotted line in Fig. 2. We can see the GRB 
rate is incompatible with the expectation from the ordinary SFR. The solid line shows the 
cumulative distribution of the 72 Swift GRBs. The result from an evolving IMF suggested 
by Dave (2010) is shown as the solid line, which agrees with the observed data very well. 



^ See [http : / /swift . gsfc . nasa.gov/docs /swift / archive / grb_table 
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From Eqs (8), (11) and (14), we can obtain the observed number of GRBs in an ob- 
server's time interval Atobs? and with Q in the interval Q — {Q + dQ), 



AiV(Q) 



^.^^ ^grb(^) AnsdVjz) dz 
LunM 1 + ^ 47r dz dQ 



dQAU,, (13) 



where Lii^{z) is determined by Eq. (fTTj) . 

We use the luminosity cut = 0.8 x lO^^erg s~^ (Li 2008). The cut in luminosity and 
redshift minimizes the selection effect in the GRB data. The total number of GRBs with 
Liso > -Z^iim is 72. The distribution of Q for the 72 Swift GRBs is plotted in Fig. 3. Because 
of the flux limit of the detector (Kistler et al. 2008), the model deviates from the data at 
z > 4. The sohd line shows the best fit of the i?GRB to the first six data points using an 
evolving IMF. The dashed line shows the best fit from non-evolving IMF. 

Fig. 4 shows the distribution of Q for all 122 Swift GRBs in the sample. The solid line 
is N{Q) calculated by Eq. [T3]with the normalization and the luminosity function parameters 
determined above, and calculated by Eq. [TTl We can see that the modeled N{Q) fits 
the observational data very well with Xr = 1-14. However, there is an obvious excess in the 
number of GRBs in the bin of < Q < 1, which might be caused by statistical fluctuations. 
If we exclude it, the Xr decreases to 0.32. 



5. Conclusions 

In this Letter, we have presented that the redshift distribution of Swift GRBs with 
measured redshifts and calculated luminosities can be successfully fitted by the SFH with 
an evolving stellar IMF. It is widely considered by current theories that only massive stars 
with masses larger than the critical value can produce long GRBs. The evolving stellar IMF 
becoming increasingly top heavy at larger z suggested by Dave (2010) can lead to more 
GRBs produced at high redshifts. 

Kistler et al. (2008) considered several possible reasons for the discrepancy between 
the Swift GRB rate and the SFH. They showed that the Kolmogorov-Smirnov test does not 
favor an interpretation as a statistical anomaly. Selection effects are also unlikely to cause 
an increased efficiency in detecting high-redshift GRBs. Although Kistler et al. (2008) have 
argued that alternative reasons are possible (e.g., evolution in the fraction of binary systems, 
an evolving IMF of stars, cosmic metallicity evolution), they did not give a quantitative 
analysis or a detailed discussion of the evolving IMF. We enlarged the GRB sample with 
122 long GRBs and used a reasonable evolving IMF. The results in this paper indicate 
that the evolving IMF may explain the redshift distribution of Swift GRBs. If the redshift 
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distribution of GRBs and SFH are well measured, GRBs would be used to probe the stellar 
IMF at high redshifts. 

We thank D. M. Wei for valuable discussions and an anonymous referee for useful sug- 
gestions and comments. This work is supported by the National Natural Science Foundation 
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Fig. 1. — Distribution of the isotropic-equivalent luminosity for 122 long-duration Swift 
GRBs. The solid line is plotted according to Eq. 12 using best fitted parameters. 
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Fig. 2. — The cumulative distribution of 72 Swift long GRBs with Ljso > 0.8 x 10 erg s 
(stepwise solid line). The dotted line shows the GRB rate inferred from the star formation 
history of Hopkins & Beacom (2006). The dashed line shows the GRB rate inferred from 
star formation history including an evolving IMF. 
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Fig. 3.— Distribution of Q for 72 Swift GRBs with Li^o > 0.8 x lO^^erg s"^ (the sohd 
histogram, with the number of GRBs in each bin indicated by a dark point with Poisson 
error bars). The solid line is the best fit of the GRB rate. The dotted curve shows the best 
fit by non-evolving IMF. 
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Fig. 4. — Distribution of Q for all the 122 Swift GRBs, the points with error bars represent 
the number of GRBs lying between ~ + 1. The sohd curve shows the best fit by Eq. 
13. The dotted curve shows the best fit by non-evolving IMF. 



